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Abstract. Quantum oracles play key roles in the studies of quantum 
computation and quantum information. But implementing quantum or- 
acles efficiently with universal quantum gates is a hard work. Moti- 
vated by genetic programming, this paper proposes a novel approach 
to evolve quantum oracles with a hybrid quantum-inspired evolutionary 
algorithm. The approach codes quantum circuits with numerical values 
and combines the cost and correctness of quantum circuits into the fit- 
ness function. To speed up the calculation of matrix multiplication in the 
evaluation of individuals, a fast algorithm of matrix multiplication with 
Kronecker product is also presented. The experiments show the validity 
and the effects of some parameters of the presented approach. And some 
characteristics of the novel approach are discussed too. 

1 Introduction 

Quantum computation is a flourishing research area and it has been believed 
that quantum computers hold a computational advantage over classical ones PQ. 
Generally, for simplification in the studies of quantum computation, a transfor- 
mation or even a quantum gate which isn't directly implemented physically, is 
treated as a black box, i.e. quantum oracle. However, one of the challenges im- 
plementing practical quantum computers is to design the quantum oracle with 
quantum circuits made of available quantum gates . Mathematically, designing 
a quantum oracle can be formulated as decomposing the expected unitary ma- 
trix to some smaller matrices which correspond to the primary quantum gates. 
But the best known mathematic algorithms are not very efficient P^. Moreover, 
the potential mechanics of quantum computation is not well understood yet. So 
it's very difficult to develop heuristic approaches designing quantum circuits at 
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present. Finally, the designed quantum circuit may be inefficient. Appropriate 
optimization techniques are required to reduce its cost, such as rewriting rules 
based |3| or templates based 0] techniques. 

However, we could get some surprised results with genetic programming 5 . 
Genetic programming is a sort of robust optimization algorithm, which mimics 
natural evolution that encodes the solution with chromosome, crossovers the old 
individuals according to the fitness, mutates them with probability and obtains 
the new individuals each generation over and over again. It is used to figure 
out complex optimization problems and requires little prior knowledge of the 
problems. All the advantages of genetic programming smooth out the above 
mentioned difficulties in designing quantum oracles, thus it has been used to 
automatic designing quantum circuits [6171819111)11 1I12| . 

In this paper we develop a novel approach to evolve quantum oracle with a 
Hybrid Quantum-inspired Evolutionary Algorithm (HQEA) which has been suc- 
cessfully applied in numerical optimization and 0-1 knapsack problems (see 
for detail introduction of HQEA). Taking the matrix corresponding to an oracle 
as input, the presented approach can achieve designing and optimizing the quan- 
tum circuits at the same time. This paper is organized as follows. In Section [21 
the novel approach evolving quantum oracles is presented. To speed up the eval- 
uation of individuals, a fast algorithm for matrix multiplication with Kronecker 
product is presented in Section [3] too. The validity of the presented approach 
is shown and the effects of some parameters of the approach are discussed in 
Section 0] In the last section, some conclusions are drawn. 

2 Novel approach evolving quantum oracles 

The novel approach takes HQEA as optimization algorithm. To take advantage 
of HQEA, it is necessary to encode quantum circuit with numerical value, which 
is different to the symbol notations in the previous works. Here, the different 
quantum gates or the same gates operating on different qubits are treated as 
different cases, which correspond to different number. Although it is possible that 
several quantum gates are operated parallclly, for simplification, it is assumed 
that only one quantum gate is applied on some qubits each time. 

In the previous works (3141617181911011 1115] , the quantum gates operating on 
non-adjacent or multiple qubits are taken as primitive gates. The commonly 
used non-adjacent quantum gates, for example, are generalized CNOT in which a 
NOT operation on a qubit is controlled by a non-adjacent one. In fact it's difficult 
to implement the quantum operations on distant qubits. It's more reasonable 
to take some one-qubit and adjacent two-qubit operations as primitive gates. 
One might argue that a realizable circuit could be obtained by replacing the 
generalized CNOT with the equivalent circuit composed by adjacent CNOT. 
But it is likely to generate redundances, see Fig. as an illustration. In this 
paper, the available primitive gates include Phase(S), 7r/8(T), Hadamard(H) and 
CNOT p. But our approach is not confined to these gates. Generally speaking, 
if ni one-qubit gates and ri2 adjacent two-qubit gates are available, there are 



N = ri\m + 2n2{m — 1) + 1 different cases on m qubits (including the quantum 
wire). So k = |~log 2 N~\ bits are required to encode one case. And if a codon is 
decoded as integer s, it corresponds to the case indexed by [%ir J • An individual 
representing a quantum circuit is consisted of gk qubits where g is the maximal 
number of allowable gates for the circuit. 




Fig. 1. Replace non-adjacent CNOT with adjacent CNOT 



Since even the same primitive gates have different costs depending on re- 
alization technologies of quantum computers, the cost of the designed quan- 
tum circuits should be considered. However, our approach is independent to 
the used cost function. In our work we just assign the one-qubit gate costs 
one, two-qubit gate costs two, and the quantum wire costs nothing. The cost 
of a circuit, allcost, is the total cost of the gates presented in the circuit. 
Let C be the implemented circuit which corresponds to the matrix A(C) = 
(/ 2 ™9 <£> A g ® I 2 "g ) x • • • x (J 2 m i <S> Ai <g) /2«i ), where m l and nt (i = 1, • • • , g) 
are respectively the numbers of qubits before and after the gate Ai. If G be the 
goal matrix the circuit should implement, the correctness of C is defined by 
such as 

|tr(G+A(C))| 

correctness(G) = - — (1) 



The fitness function takes allcost into account as well as the correctness of the 
circuit. Evolving quantum oracles is time-consuming when the scale of problem is 
very big. Sometimes we could be satisfied with a non-optimal circuit whose cost 
is less than some value. Such a value is called the satisfying cost (satcost). But 
it's more important to obtain a correct quantum circuit. So a tradeoff between 
the cost and correctness should be taken. We take two thresholds award and 
punish, and then the fitness function is defined as 



fitness(C) = award x (allcost — satcost) +punish x (1 — correctness(G)) (2) 



As an evolutionary algorithm, the termination condition of the proposed 
approach is meeting the satisfying cost as well as fulfilling the correctness, or 
evolving allowable generation. 



3 Fast matrix multiplication with Kronecker product 



Lots of matrix multiplications are required in the process of evaluating indi- 
viduals. It is well known that we should perform 0(n 3 ) multiplications when 
two n x n matrices are multiplied naively. Although the best algorithm currently 
known has an asymptotic complexity of 0(n 2 376 ) |14j . some improvement is still 
possible for the matrix multiplication with Kronecker product which is required 
in this paper. Firstly, we put up some conventions. 

— Kronecker product, matrix multiplication and scalar product are denoted as 
(g>, x and • respectively. The priority of • is higher than others and x takes 
on the lowest priority. 

— A n k ^ denotes arbitrary n x n matrix A treated as k x k blocks, each of which 
is a j x ^ matrix. Specially, A n means An and l n denotes the nxn identity 
matrix. 

— B(i,j) n k denotes the block located at the i-th row and the j-th column in 
^mlk- I n Particular, A(i,j)i denotes an element of A, for short, a^-. 

The trick is based on the block multiplication and the calculation is shown 
as follows. 
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D (p,<l)k = ^2 {a p i ■ B(i,j)nk(l,q)k) (p,q=l 



1=1 



(5) 



It's very clear that obtaining D(p,q) requires 0(nk 2 ) multiplications, and 
then 0(n 3 k 2 ) multiplications to obtain D and lastly 0(m 2 n 3 k 2 ) multiplications 

to obtain F. So the fast algorithm speeds up ^™^ 3fc2 = (mfc)°- 376 /n - 624 to the 



original best algorithm. Simple algebra shows that our fast algorithm exceeds 
the best traditional algorithm when the number of qubits before the gate, M = 
log 2 m, the number of qubits after the gate, K — log 2 k, and the number of 
qubits the gate processes, N = log 2 n, satisfy M + K > 1.66-/V. 



4 Experiments and discussions 

The presented approach is not confined to circuit design or circuit optimization, 
and the key factor is the satisfying cost. When this constraint is loose correspond- 
ing to a bigger satisfying cost, the algorithm performs mainly as an automatic 
designer which aims to discover a circuit implementing the desirable unitary 
transformation, and regardless whether the result is optimal. If this constraint is 
tight corresponding to a smaller satisfying cost, the algorithm not only tries to 
dig out a circuit functioning as desired, but also to reduce the cost of the circuit. 
As an instance, the optimal circuit of oracle SWAP is CNOT * CNOT2 * CNOT 
(CNOT2 denotes the CNOT taking the latter qubit as control bit and the former 



qubit as controlled bit) such as Fig. 2(a) which costs 6, but if we assign satcost 



as 8, another circuit could be obtained by our approach, such as Fig. 2(b) 
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(a) Optimal circuit 



(b) Non-optimal circuit 



Fig. 2. Equivalent circuits of SWAP 



Table 1. Comparison with various satisfying costs 



Goal oracle Optimal 


Satisfying 


Maximal number Maximal 


AS 


ST OT a 


cost 


cost 


of allowable gates 


generation 






entangle2 3 


4 


6 


100 


86.1 


20 4 




6 


6 


100 


14.8 


20 


entangleS 5 


6 


8 


200 


141.7 


20 10 




8 


8 


200 


48.65 


20 1 


controlled-B 7 


8 


8 


500 


111.5 


20 3 




10 


8 


500 


62.5 


20 1 



" AS = average generation of success, ST — times of success in all tests, OT = times 
of getting optimal results in all tests. 

To study the affection of satisfying cost on evolving results, we apply novel 
algorithm on the oracle entangling two qubits, denoted by entangle2 as Fig. 



the oracle entangling three qubits, denoted by entangle3 as Fig.^Jand the oracle 
implementing controlled-phase, denoted by controlled-S as Fig. [3] Each case is 
tested 20 times and in all the test, the number of quantum chromosome is 20, 
measurement times of each quantum chromosome is 10. The comparison results 
are shown in Tabled In the experiment, all the tests obtain the circuits satisfying 
condition successfully. However, obtaining the circuits meeting more rigorous 
constraints needs more evolving time. Additionally, it can be found that there 
are many equivalent circuits implementing the same oracle, although some of 
them cost differently. Some of these circuits are illustrated in Fig. I3I4I5I 



(a) Optimal one(cost 3) 



(b) One of non-optimal circuits(cost 
5) 



Fig. 3. Circuits entangling two qubits 




(a) Optimal (b) One of non-optimal cir- 

one(cost 5) cuits(cost 7) 



Fig. 4. Circuits entangling three qubits 



(a) Optimal one(cost 7) 



(b) One of non- 
optimal circuits(cost 
8) 



(c) One of non- 
optimal circuits(cost 
9) 



Fig. 5. Circuits implementing controlled-phase 
To effectively discover desired quantum circuits for different cases, it is useful 
to adopt the appropriate reward-punish factor in the fitness function. Table [5] 
shows the evolutionary results of oracle entangle2 with different reward-punish 
factors, where each case is tested 20 times. It is found that bigger punishment 



to the error of circuits is required to get the correct circuits, with the same 
satcost, g and other parameters. Another fact is that with the same g and 
other parameters, to get the correct circuits bigger punish is required for bigger 
satcost. 

Table 2. Comparison with different reward-punish factors (reward = 1) 



Satisfying 


Maximal number Maximal 


ST 


AS" 


punish 


cost 


of allowable gates generation 
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6 100 








1 




100 


20 


16.9 


5 




100 


20 


10.65 


20 




100 


20 


12.8 


100 




100 


20 


15.55 


1000 
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8 200 








1 




200 








5 




200 


20 


32 


20 




200 


20 


94 


100 




200 


20 


62.65 


1000 


10 


8 500 








1 




500 








5 




500 








20 




500 


20 


40.45 


100 




500 


20 


31.1 


1000 



a ST = times of success in all tests, AS — average generation of success. 

When comparing our work with others 3 4 6 7 8 9 1Q I11I12| . we observe the 
following aspects: 

1. Problems considered: Our approach evolves quantum circuits taking the de- 
sired unitary matrix as input, while some other works are based on the de- 
scription of an oracle and evolve quantum circuits by comparing the outputs 
of the quantum oracle on random inputs with desirable ones [5]. In addi- 
tion, the simplification of known quantum circuits are considered by |3I4) . 
Notably, only reversible quantum oracles are considered in our work. Thus 
measurements, not like the cases in |7I8) . are not permitted. 

2. Primary gates set: As has been stated, only two-qubit gates on adjacent 
qubits are available in our work while non-adjacent two-qubit or even multi- 
qubit gates are taken as primary gates in other works, although it can sim- 
plify the problem of quantum circuit design. Of course, our approach is not 
confined to any special quantum gates set. 

3. Circuit representations: While all the previous works use symbolic represen- 
tations, our approach encodes the quantum circuits with numerical values. 
In despite of the apparent difference within them, all of them are equivalent. 
But it is more natural to apply evolutionary operators on the individuals 
coded by our means. 



4. Fitness function: The cost of designed quantum circuits is ignored in |6I7I8I9| . 
Moreover, the cost and correctness are dividually considered in To re- 
flect the fitness of evolved quantum circuits more accurately and expediently, 
our approach combines them together by simple reward-punish factors. 

5. Application: Our approach can be adapted to both circuit designing and 
circuit optimization. When applied to the later, more rigorous conditions 
should be assigned, such as the smaller satisfying cost and maximal gates 
number. 

6. Common challenge: The bottleneck in designing quantum circuits with evo- 
lutionary algorithms is the individual evaluation, i.e. the matrix multiplica- 
tions in the computing of fitness which have potentially exponential space 
increment and speed slow down. The intractable problem results from the 
argument that the quantum system can not be effectively classically simu- 
lated. 

5 Conclusions 

Genetic programming appears to be useful in designing quantum circuits. We 
propose how to evolve quantum oracle with a hybrid quantum-inspired evolu- 
tionary algorithm. With our approach no additional knowledge is required to 
design an optimized quantum oracle as expected. Especially, we design a novel 
approach to represent quantum circuits with numerical values in the evolutionary 
algorithm. A faster algorithm for matrix multiplication with Kronecker product 
is presented too. It speeds up the evaluation of individuals very much. Obvi- 
ously, the numerical representation and fast algorithm of matrix multiplication 
are not unique to our approach, but adaptable to other evolutionary quantum 
programming algorithms or hierarchical approaches. By assigning different pa- 
rameters, the novel approach could be inclined to designing a quantum circuit 
or optimizing it. Our approach provides insights into quickly evolving quantum 
oracles. A possible improvement to the approach may be encoding the quantum 
circuits with variable length. 
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